Analysis and Comparison of Port-Hamiltonian Formulations for Field 
Theories - demonstrated by means of the Mindlin plate 
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Abstract — This paper focuses on the port-Hamiltonian for- 
mulation of systems described by partial differential equations. 
Based on a variational principle we derive the equations of 
motion as well as the boundary conditions in the well-known 
Lagrangian framework. Then it is of interest to reformulate 
the equations of motion in a port-Hamiltonian setting, where 
we compare the approach based on Stokes-Dirac structures to 
a Hamiltonian setting that makes use of the involved bundle 
structure similar to the one on which the variational approach is 
based. We will use the Mindlin plate, a distributed parameter 
system with spatial domain of dimension two, as a running 
example. 

I. Introduction 

Distributed parameter systems described by partial differ- 
ential equations arise in systems theory from a modeling 
and a control theoretic point of view and are without doubt a 
challenging research problem, where lot of progress has been 
achieved in the last years. Also the port-Hamiltonian setting, 
originally developed in the finite dimensional scenario has 
been transfered to infinite-dimensional systems, where e.g. 
the well-known approach based on (Stokes-)Dirac structures 
(also known from the lumped parameter scenario) is avail- 
able, see e.g. [1], [2], [3], [4], [5], [6] and references therein. 

Also in mathematical physics, systems described by partial 
differential equations (pdes) are interpreted in a Hamiltonian 
setting, e.g. in [7], [8], [9] and references therein, but in 
most cases systems with trivial boundary conditions are 
considered, which is not the case in many engineerings 
applications. Therefore, the approach based on Stokes-Dirac 
structures has been setup to overcome the problem of non- 
zero energy flow through the boundary. 

A different port-Hamiltonian approach is based on a 
bundle structure with respect to independent and dependent 
coordinates (not necessarily relying on an underlying Stokes- 
Dirac structure), see also [10], [11], [12], [13], [14], [15], 
[16] which are all based on [9] but adapted to control 
purposes, i.e. modified in a sense, such that non-zero energy 
flow through the boundary can be considered (such that 
boundary ports are included) and furthermore control inputs 
on the domain and/or the boundary can be included. 

The main difference of the approach relying on Stokes- 
Dirac structures and the approach using bundles, is that the 
Stokes-Dirac scenario is based on the choice of proper energy 
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variables (flows and efforts) for which the power balance 
is formulated, whereas the second approach is based on a 
given Hamiltonian density (the total energy density) and the 
evaluation of the power balance is performed based on the 
underlying bundle formalism in order to restructure the pdes 
such that the energy flows are linked to the physics. This 
will have the consequence, that the variational derivative is 
interpreted differently and the choice of state variables is 
different, in the two mentioned approaches. 

The purpose of this paper is, that based on the weU- 
known Lagrangian setting for first order field theories i) 
the partial differential equations and the boundary conditions 
derived using a variational principle are reinterpreted using 
two different port-Hamiltonian settings, which describe the 
same physical phenomenon but using a completely different 
port-Hamiltonian representation, ii) by using the example 
of a Mindlin plate all these concepts are visualized and 
compared in great detail. 

II. Notation 

We will use differential geometric methods for our con- 
siderations and the notation is similar to the one in [7], 
where the interested reader can find much more details about 
this geometric machinery. To keep the formulas short and 
readable we will use tensor notation and especially Einstein's 
convention on sums. 

We use the standard symbol A for the exterior product 
(wedge product), d is the exterior derivative, J the natural 
contraction between tensor fields. By 9^ are meant the 
partial derivatives with respect to coordinates with the indices 

and [m"^] corresponds to the matrix representation of the 
(second-order) tensor m with components m"^. E.g. taking 
a second-order tensor m and a co-vector ui, the components 
of the contraction m\uj read in local coordinates as m°'^uja, 
where the summation over a is performed (Einstein conven- 
tion on sums). 

Furthermore C°°(-) denotes the set of the smooth func- 
tions on the corresponding manifold. Moreover we will not 
indicate the range of the used indices when they are clear 
from the context. Additionally, pull backs and pull back 
bundles are only stated when necessary, when they follow 
from the context they are not indicated to avoid exaggerated 
notation. 

Let us consider the bundle y ^ V, {X^,y°') -> (X^). 
The first jet manifold J^{y) possesses the coordinates 
{X^, y°',y'X), where the capital Latin indices A, B are used 
for the base manifold T) (independent coordinates) and 
denote derivative coordinates of first order (derivatives of the 



dependent coordinates with respect to the independent ones) 
as well as 

d 



dA 



A 



The jet structure also induces the so-called total derivative 

dA = dA + yldo, + y%d!^ 

acting on elements including first order derivatives and y'^^ 
correspond to derivative coordinates of second order living 
in J^'^iy), the second jet manifold. Based on the bundle 
structure y ^ V let us introduce the vertical tangent bundle 
V{y), as well as 

Ai{y)^r*{y)A{Ar*m, 

see also [7], with a typical element lu — uJady" A dV 
for Af{y) where dV^ denotes the volume element on the 
manifold V, i.e. dV^ = dX^ A... AdX'^ with dim(2?) = d 
and the functions Wa may depend on derivative coordinates. 
Furthermore, a typical element for V(3^) reads as v — v°'da 
and when depends on derivative coordinates we call v a 
generahzed vertical vector field, see [9]. 

III. Background material 
A. Geometric preliminaries 

We will consider densities ^ in the sequel (a quantity that 
can be integrated), J = FdV with J" 6 C°°(ji(A')) (we 
restrict ourselves to the first-order case). By F = J^^^ we 
denote the integrated quantity, where of course a section of 
the bundle y V, i.e. a map y — leading to yA — 

dA^{X) has to be plugged in to be able to evaluate the 
integral properly. 

Proposotion 1: Given the density ^ = J'dV and a gen- 
eralized vertical vector field v : V ^ together with 
its first jet-prolongation — v°'da+ dA{v°')d^, see [9], 
[7], we obtain the decomposition 



(1) 



IV Jv JdV 

Here the map 5^ ~ SaJ'dy'^ A dV^ (corresponding to the 
Euler Lagrange operator of see [7], with the coefficients 
6aJ^ = daJ- — dAd^J- (called the variational derivatives) is 
used, as well as the boundary operator 5^5 — d^J-dy°' A 
dVA with dVA = dA\dV (the boundary volume form). 

Proof: The proof follows by evaluating the Lie- 
derivative of the geometric object J with respect to the vector 
field j^{v) 

j\v){TdV) = {v''{d^T-dAd^T) + dA{v''d^T))dV 



{v'^6^F + dA{v"d^F)) dV 



(2) 



and applying the Theorem of Stokes [9] to (|2]i. ■ 

The relation ([T]l will be of key interest in the forthcoming, 
since it provides a natural decomposition of the expression 
/u i^{'v){FdV) into a term on the domain T) and one on the 
boundary dV. Important is the case when the generalized 
vector-field v is Unked to the solution of a pde system (via 



its semi-group, that v may generate), then the formal change 
of F = FdV along solutions of a pde system can be 
computed as i^{v){FdV) (provided all operations are 
admissible), which we denote by F = i^{v){FdV) in 
this special case. 

B. Dirac structures 

Based on the space of power variables T x £ (flows and 
efforts) and the symmetric bilinear pairing 



< (/i,ei),(/2,e2) >:= (ei,^) + (e2,/i) 



(3) 



where (•, •) is the dual product of the linear spaces and 
£ = F* a Dirac structure is a linear subspace D C Tx£ such 
that D = with respect to the pairing ©. For (/, e) € D 
one has (e, f) — such that the Dirac structure preserves 
power This concept can be transfered to the case where 
and £ are spaces of vector-valued functions over a spatial 
domain V, then infinite dimensional systems are the focus, 
and to allow for non-zero energy flow through the boundary 
the so-called Stokes-Dirac structure is introduced, see [1], 
[2], [3], [6] and section EH 

IV. Lagrangian Framework 

In this section we recapitulate the well-known Lagrangian 
framework for first-order field theories, and we will derive 
the partial differential equations as well as the boundary 
conditions in a geometric fashion. Thus, we consider a 
bundle 

Q^Vl, {q'',t\X^)^{t\X^) 

where we use the shortcut x' = such that the 

independent variables are the time t° and the spatial ones 
X^ , A = 1, . . . dim(2?i) — 1 = n . A first order Lagrangian 
takes the form 



£: J'(Q)^"aV*(Pl) 



(4) 



£ = £w with C e C°° [J^ (Q)) together with the volume 
element uj that meets 

w = dt° A dX^ ... A dX" , = d^\uj , i = 0,...,n 

The variational problem for a section s : T>l^ Q is the 
followindl 



= 0, 



(5) 



e=0 



where the flow V'e is used to deform sections s : Vl Q and 
whose generator is a vertical vector field vl : Q ^ ^(Q). 

This is a well-known problem and treated for example 
in [7], [9] and references therein. It is obvious that Q is 
equivalent to 



(6) 



'At this point tlie pull back is essential and therefore indicated, i.e. 
(ipcos)* S, means, that the first prolongation of (tpios) has to be plugged 
in into £ = Cu! in order to evaluate the integral properly. 



see for example [7] and based on ([T]l where we replace T 
be L we obtain the decomposition 

/ i;lJ(5£+ / i;lJ^^£ = 0. (7) 

Consequently, the partial differential equations for a first 
order Lagrangian follow as 

d^{c) = o, (8) 

and the boundary term is the second term in ^ and reads 
in local coordinates as 

/ vldlCuj, = 0. (9) 

JdVL 

The boundary conditions can be fulfilled by either allowing 
for no variations on (a part of) dV^, i.e. v£ = or by 
(9^£ = or by a combinations of both approaches. Possible 
is also the inclusion of external boundary variables ^ such 
that 

/ v2{dl,C-FUuj,^0. (10) 

has to be met. 

V. Port-Hamiltonian Picture 

Now we turn to the Hamiltonian picture, where we dis- 
cuss two different port-Hamiltonian formulations. We will 
restrict ourselves to systems without dissipation and without 
distributed control for simplicity, but these properties can be 
included in both formalisms in a straightforward manner, see 
[2], [11]. 

A. Geometric approach based on underlying bundle struc- 
ture 

We will introduce port-Hamiltonian systems described by 
pdes based on a power balance relation, such that the power 
balance relation together with the structure of the equations 
represent the physical process. 

Definition 1: A port-Hamiltonian boundary control sys- 
tem without dissipation on a bundle X ^Vn, ^ 
[X"^) takes the form of 

X = J {5^) (11) 

with the Hamiltonian ^ UVl , D. ^ dX^ A ... A dX'^, 
where Ti. G C°° {J^ [X]) and additional boundary conditions 
(possibly including boundary inputs, optionally leading to 
so-called boundary ports). The map J is of the form J : 
h.f{X) V{X) where ^7 is a skew-symmetric map. 
In general the map J can be a differential operator, see our 
paper [16], but within this contribution we exclude this case 
(since in many examples, e.g. mechanics this is not required). 

Now we make use of proposition [T] and replace by in 
([U. Setting w = i we obtain 

H ^ I x\5^^^ [ x^d^nilA (12) 

where Ha = Oa] H, which reflects the power balance, since 
the total change of the functional H along solutions of STT[ . 
is affected by a boundary port (if it exists) depending on the 
boundary conditions. See e.g. [11] for a formal introduction 
concerning the boundary ports. 



B. Approach based on underlying Stokes-Dirac Structure 

Following [6] we shortly recapitulate the port-Hamiltonian 
framework based on Stokes-Dirac structures. For more de- 
tails we refer to [1], [2]. 

We consider the space of flows J' and the space of efforts 
£, which are spaces of vector-valued functions over a 
spatial domain V. Given J'sd a skew-adjoint matrix 
differential operator [6], the space 

B = {if,e,w)eTx£ xW\f = -JsDe,w ^ Bv{e)} 

is a Stokes-Dirac structure, regarding the pairing 

< (/i,ei,u;i) , (/2, 62,-102) >= 

= JT,[e1f2 + elfi\dV + Bjsr,{wi,W2)dA 

where B j^^ is a boundary differential operator induced by 
JsD any by slight abuse of notation corresponds to the 
boundary volume element. The map B-d is a boundary 
operator and the boundary variables are w. 

If (/, e, w) G B then 

= f e'^fdV+l f BjsA'^,w)dA (13) 

holds. Given an energy density T-LdV where % depends on 
the energy variables, a port-Hamiltonian boundary control 
system without dissipation can be stated as 

f = -JsDe, w^Bv{e) (14) 

where the energy variables /, are linked to the state variables 
X via / — ~x and the efforts variables follow from e = 

Remark 1: Originally, in [6] instead of e — d^H the 
authors use e — 6^H where H — Jjy T-LdV, but since 
TL depends on energy variables, the variational derivative 
degenerates to a 'partial' one. 

Furthermore, from H ~ HdV and the relations (ITSl l and 
([l4l i one has 

H^l f Bj,^{w,w)dA (15) 

since H ^ d^UidV = - e^fdV. 

Remark 2: To derive this energy balance also proposition 
[T] can be applied, but since no jet-variables are included it 
simplifies to H = J^j^{v){ndV) = J^v°'da'HdV and to 
derive (T5[ a further integration by parts must be performed, 
since v corresponds to — / which involves the differential 
operator Jsd- 

VI. The Mindlin Plate 

Let us consider a rectangular plate with lengths lx,ly, 
where h will denote the thickness, which will be modeled 
based on the hypothesis stated by Mindlin. Therefore, we 
choose as independent coordinates the vertical deflection w 
of the mid-plane as well as the rotations of a transverse nor- 
mal to the X and Y direction termed 1/; and (p, respectively. 
The kinetic energy density /C and the potential energy density 
V can be stated as 



/c 



and 



V 



-kGh [{wx - -0)^ + {wY 



2 2 

+ \iD{i^X + V(t)YiJx) + £'(0y + vMx)), 

2 

where v is the Poisson ratio, k = and G, I? are the 
plate stiffness and the plate module, respectively, see [6] and 
references therein. 

Remark 3: The subscripts t, X, Y correspond to the 
derivatives with respect to these independent variables, ac- 
cording to the jet-bundle structure in the Lagrangian frame- 
work. The subscripts x, y to be used later, correspond to 
quantities which are connected to the spatial variables X and 
Y but they must not be confused with derivative variables. 
Since we are in a time-invariant setting, we will use later 
on also the ' notation, for time derivatives, instead of the 
subscripts t. 

To derive the equations of motion we will use the variational 
principle in a Lagrangian setting. Then given the partial dif- 
ferential equations, we will interpret them in a Hamiltonian 
setting, either using the approach presented in section IV-AI 
and using an approach based on the Stokes-Dirac structure 
as in section I^B] 

A. The Lagrangian picture 

In the Lagrangian framework we consider the bundle 

Q^VcAw.i^,<P.t.X,Y)^{t,X,Y) (16) 

together with the Lagrangian density £ = Cuj with C = 
/C — V, uj — dt AdX AdY. The variational derivatives follow 
form the chosen bundle structure ( fT6b and follow to 

Sw = dm- dtdl - dxd^ - dYd^, 



w 



w 



S-ip — — dfd^ — dxd^ — dYd^ 
6^ = ~ dtdl - dxdf - dYdJ . 

From Si^C ~ 0, Sjf,C = 0, J^^C — corresponding to (O we 
derive the partial differential equations 

phwtt = kGh{wxx - i/'x) + kGh{wYY - 4>y) 
Pini'tt = kGh{wx -4') + D{il;xx + v4ixy) 



12 



+ ]^D{1 - v){i;YY + (t>XY) 



PT7:4>tt = kGh{wY - (t>) + D{(j)YY + i^Tpxy) 



12 



D{iPx + v<Py) = -d^C 



If we introduce 

Af, = 

My = D{^Y + v^x) = -dlC 

M^y = D^^-^{4jY + 4>x) = -dlC = -dfC(n) 

= kGh{wx~i^)^-d^C 
Qy — kGh{wY ~ (f>) — —d^C 

then the equations of motion take the familiar forrrQ 

phwtt = dxQx + dYQy 
Pj^At = Qx + OxM^ + dyM^y (18) 

Pj^'Ptt = Qy + dYMy + dxM^y. 
The boundary conditions follow from 

{wLdlC + i^Ld^C + cj)LdlC) d,\ [dt AdX AdY) = 

dVc 

where the variational vector field vl takes the form vl = 
WLdw+ipLd^ + ipLd^ and on the time-boundary no variation 
takes place (i.e. when i = then vl = 0) and (VT\ has to 
be used. Therefore we have 



+ ^lM^ + (j)LM^y) dt A dY 
iwLQy + ipLM^y + (j)LMy)dt AdX ^0 (19) 



such that, e.g. if at X = we have that wl is arbitrary, then 
Qx has to vanish or has to be compensated by an external 
boundary term as in ( fTOb . such that the familiar boundary 
conditions are recovered. 

B. The Hamiltonian picture 

Based on the partial differential equations ( fTsT i and the 
boundary conditions ( fT9b we discuss the two presented port- 
Hamiltonian formulations as well as the power balance 
relations corresponding to the particular representation. 

1) Geometric approach: Now we consider the bundle 
(which is different form the one in the Lagrangian setting) 

X -> T^H, {■w,ip,(f>,Pw,P4„P<t>,X,Y) (X,Y). 

From the Legendre transform we derive the temporal mo- 
menta ^ d^jC, p^ = d^pC, p^ = d^C which read as 

Pni = phwt , p.0 = PTTrV't , P4> = P— (20) 



12 

and the Hamiltonian follows as 

T-L = wpw + ipp-4, + 



12 



L. 



+ -D{l-y){iPxY +<^>xx)■ 



\v^ the coordinates {w,ip, (f>,Pw,Pt(>,P4,) together with (|20] | 
one has -H = /C + V. 

-In order to be comparable with the hterature we sometimes use dxaivi 
dy although in a strict mathematical sense it should be dx and dy ■ 



To derive the Hamiltonian formulation as in (fTTT i we set 

X = {w,il>,<l),Pw,Pi),P4,) and obtain 

which reads as 



(21) 



w 
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The variational derivatives in this setting take the form 



dxd^ 
dxd^ 



-dydl 
dydl 



and 5p^ = dp^ , 5p^ = dp^ , 5p^ = because of the differ- 
ent bundle structure compared to the Lagrangian approach. 
The boundary ports follow from ( fT2] ) and we obtain 



H 



x"d^ndA\ [dX A dY) 



dVi-_ 



{x'^d^HdY - x'^d^HdX). 



From the special choice of the Hamiltonian H we observe 
that the expressions d^C and d^H correspond (apart form 
the sign), and therefore we also have 





— Qx 




— Qy 








— Mxy 











which consequently leads to the power balance relation 



H 



dQx + i>M.x + ^Mxy)dY 
{wQy + TpMyx + ipMy)dX (22) 



which is based on (fT2l i. 

Remark 4: It should be noted that depending on the 
boundary conditions, along dV a boundary port appears only 
if in the pairings x^d^H both 'players' are not equal to zero. 
Furthermore it should be noted that in i" only 
remain, since in H only jet variables with respect to w, ip, (j) 
appear, i.e. there is a wx present but no {pw)x and so on. 

2) The Stokes-Dimc approach: This approach is not 
based on a bundle structure, which distinguishes dependent 
and independent variables strictly, but uses so-called energy 
variables instead. Therefore, let us introduce the strain vari- 
ables as [6] 



-i^x 

~4>y 

-{ipy + cj)x) 
Wx — ^ 
Wy — (j). 



(23) 



Then one can introduce as state x which consists of the 
momentum variables and the strains 



X = {phw,T, 
and the Hamiltonian % 

P 



r ^1 ^ 

xz, 1^ yzi P ^^Vt P 



; 1^X5 Ty , xy) 

/C + V can be rewritten as 

+ hw^) 



and 

V = 
such that e 



2 {Qx^x 



Qy^y^ ^^xy^ xy ^^x^x '^-^y^y) 



d^H follows to 



(w, 



}y,lJj,4),-Mx,-My,-Mxy). 



The partial differential equations can be stated as 

X = Jsoe = Jsod^T-L 
where Jsd takes the form 



(24) 






dx 


dy 

















dx 








-1 














dy 











-1 














1 











-dx 





-dy 








1 











-dy 


-dx 











-dx 


























-dy 




















-^dy 


-dx 












The energy balance follows from equation (flSl l where (fTsT i 
has to be evaluated by a further integration by parts, since 
Jsd is a differential operator which contributes to the 
boundary expression, i.e. Bj^^ has to be constructed, see 
[6]. 

VII. Discussion and Comparison 

In this section we will discuss the main differences of 
the two presented port-Hamiltonian scenarios as in sections 
(IV-Ab and (IV-Bb where we will highlight these aspects by 
focusing on the presented example, the Mindlin plate. 

A. State variables 

The state variables in the geometric approach consist of 
the displacements/deflections and the temporal momenta, i.e. 

Pw ; Pip 1 

Pff,) in our example where in many 
cases (mechanical systems) the temporal momenta can be 
derived from a given Lagrangian by means of the Legen- 
dre transformation. These temporal momenta are introduced 
mainly to obtain explicit partial differential equations where 
the state variables are differentiated with respect to a curve 
parameter, which is the time. 

In contrast to this, in the approach based on the Stokes- 
Dirac structures, energy variables are used, such that the 
strains are introduced in mechanical applications and one 
has X = {phw,Txz-,Yyz,Pj2^,Pj^(l>,Yx,^y,Yxy)- From a 
conceptional point of view the use of energy/power variables 
may be beneficial since they are linked to the power balance 
relation in a simple manner, but as can be seen already using 
the MindUn plate example, the five strain variables have to be 



derived from the three independent deflection/displacement 
variables {w, V', 0) by differentiation, see (l23l l. such that 
additionally to the partial differential equations (|24] | also 
the compatibility conditions ( l23T l must be listed, such that 
a constrained Hamiltonian representation is apparent. 

B. Control issues 

From a control point of view control methods like damping 
injection or control by interconnection can be performed 
equivalently using the two presented Hamiltonian represen- 
tations, see for example [14], [5] where the Timoshenko 
beam is analyzed. However as stated also above, the use 
of energy variables allows for controlling the system for 
instance to zero strain configuration, but the global position 
in space cannot be controlled in a straightforward manner 
since the deflection/displacement coordinates do not enter 
the formalism, in contrast to the approach as in ( IV-Ab . 

C. The skew-symmetric operators J and Jsd 

By inspection it becomes apparent that J and Jsd differ 
significantly, since J is no differential operator in contrast to 
Jsd- This also has severe consequences for the expressions 
as in (fTTT i and e = as in (fT4l i. where it is vice 

versa, i.e. (5 is a variational derivative and in the Stokes- 
Dirac approach a partial derivative appears. Let us consider 
for instance the fourth equation of ( l24l l which reads as 

p—^ = 1 • - dx{.-M..) - dvi-M^y) 
as well as the fifth equation of (1211 1 which is 

Pi, = SibU = -{d^ - dxd^ - dYdJ,)H. 

From d^H = —Qx as well as from 9^7^ = , d^H — 
Mxy we easily observe that a part of the variational derivate 
— — dxd^ — dyd^ is incorporated in Jsd (namely 
(1, ~dx, —dy)) and that the variables Qx, Mx and Mxy are 
used directly. 

D. The energy balances 

From the relation (fT2l l the power balance is derived easily 
once the state x as well as the Hamiltonian density Ti. is 
chosen - this is very simple, since J is no differential 
operator and does not contribute to the boundary term. This 
is different in the approach as in section ( IV-BI ) since the 
boundary operator Bj^^ has to be derived from the special 
choice of Jsd and by an additional integration by parts one 
ends up again by the same relation as in (l22l i. 

VIII. Conclusion 

We have presented two different port-Hamiltonian rep- 
resentations based on a given set of partial differential 
equations together with their boundary conditions derived 
by the Lagrangian formalism using jet-bundles. By means of 
the running example, the Mindlin plate, we have extensively 
discussed and compared these two different Hamiltonian 
scenarios. Further investigations should also include the field 
theoretic Hamiltonian concepts coming from mathematical 
physics, like the polysymplectic and/or the multisymplectic 
approach as in [7], [8] or in the spirit as in [17]. 
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